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Solution of Linear Systems by a Singular Perturbation Technique 


MARK D. ARDEMA 
Ames Research Center 

Abstract — An approximate solution is obtained for a singularly perturbed 
system of initial-valued, time- invariant, linear differential equations with 
multiple boundary layers. Conditions are stated under which the approximate 
solution converges uniformly to the exact solution as the perturbation param- 
eter tends to zero. The solution is obtained by the method of matched asymp- 
totic expansions. Use of the results for obtaining approximate solutions of 
general linear systems is discussed. An example is considered to illustrate 
the method and it is shown that the formulas derived give a readily computed 
uniform approximation. 


I . INTRODUCTION 

Linear systems of ordinary differential equations frequently occur in 
applications of control theory. Although the general solutions of such 
systems are well known, the computation of solutions for specific systems can 
be laborious if the order of the system is high. In applications in which it 
is desired to keep computational effort low and in which extreme accuracy is 
not required, it is appropriate to seek approximate solutions. The approxima- 
tions that result in the greatest simplification are those in which derivative 
terms are neglected. Singular perturbation theory provides a systematic way 
of studying such approximations and coping with their shortcomings. 
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Before proceeding, some notation will be introduced. Given an n * n 


matrix A with elements and an n- dimensional vector x with components 

x ^ we define 

*a£ - matrix formed by taking i— through j— rows of A and k^ 
til 

through l— columns 

jX ■ vector formed by taking i— through j— components of x 

We note that *A^ ■ a,.. ^A* = A, oc = x J , and *x - x. The determinant of 

i j ij n n i i n 

i k i k 

_jA^ will be denoted The following identity is true for i = 1, •••, 

n - 1 and will be needed subsequently: 


V - i+1 D i+1 
n n n n 


i. i+1 1 

o _ A 

f i+1 A i+1 l 

|-1 

i+M] 

“ ♦ , A 

xi in 1 

l nn j 


n ij 


( 1 . 1 ) 


We are now ready to formulate our problem. Consider the order system 


Hx 

E — = Ax; x(0) = x 
at o 


( 1 . 2 ) 


where A is a constant matrix and E is a matrix with components e^ such 
that 


e ij * 0; 1 * 1 


e li - 1; i - 1, 2, 


m 


(1.3) 


'ii 


i-m 

e ; i - m + 1 , 


n 


where e is a small parameter. Although the system (1.2) has a rather special 
form, its solution may be used to obtain approximations to more general systems, 
as discussed in section IV. Assuming n to be large, it is our goal to obtain 
an approximate solution to (1.2). 

Since e is small, it is natural to seek an approximate solution of 
(1.2) by setting e * 0. The resulting system is 


r; njoPueimuTY of tip. 
original page is poor 
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( 1 . 4 ) 


,1 ° 

m X 1*1 1 o . l.m+1 nrt-1 o 1 1 

—77 — - A x + A x ; x( 0 ) * x 

dt m m m mn n m mo 


nrt-lg nrt-1. 1 1 o nrt-1 m+1 m+1 o 

* ® A X » A X 


n 


n m m 


n n 


n 


and is called the reduced system. We call the components of ^x the "slow" 

jjj+1 

variables and those of x the "fast" variables. It is obvious that, in 

n 

general, the solution x°(t) of the redi-ced system will not satisfy all of 

the initial conditions and thus, at least locally, the behavior of the solution 

of the reduced system will be radically different from the solution of the 

exact system x(t). For this reason, we call ( 1 . 2 ) a singularly perturbed 

system. The best that can be hoped for is that ^x° is a good approximation 

of *x uniformly in the domain of interest and that ^^x 0 is a good ap- 
m n 

proximation of ^^x everywhere except near t = 0 . The conditions under 
which this is true will be established subsequently. 

Singular perturbation theory is concerned with the relation between the 
solutions of the exact and reduced systems of singular perturbation problems 
and with constructing asymptotic series representations of the solution of 
the exact system. The fundamental results of the theory are due to Tihonov 
[l] and Vasileva [ 2 ]. Reviews of singular perturbation theory and methods 
are found in the books of Wasow [ 3 ] and O'Malley [ 4 ]. Singularly perturbed 
linear systems have been studied by several authors [c.f. 5 ]. 

In the present treatment, solutions are developed by application of the 
method of matched asymptotic expansions. This method has been developed to 
solve certain fluid mechanics problems involving nonlinear partial differen- 
tial equations. Expositions of the method are given by Cole [6], Van Dyke 
[ 7 ], Nayfeh [8], Eckhaus [ 9 ], and Ardema [ 1 C]. In the MAE method, boundary- 
layer systems are introduced to account for the nonuniform approximation in 
the fast variables. Due to the special nature of E, there will be n - m 
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boundary layers for the system (1.2). In effect, each boundary layer, except 
the n - nr^, will be a reduced solution for the next boundary layer. Such 
systems of multiple boundary layers have been considered by Calise til, 12]. 
The reduced and boundary- layer systems are solved by asymptotic expansion 
methods. The solutions of these systems, which must be asymptotically stable, 
are then "matched" with the solutions of the reduced problem. Finally, a 
composite solution which gives uniform approximation in all variables is 
formed. In the current application, we are interested only in the zeroth 
order terms of the expansions and for this case the MAE method assumes a 
particularly simple form. 


II. THE MAIN RESULT 


Theorem 

Consider the system (1.2) and suppose that 

(-l) 1+n i D i < 0; i = m + 1, •••, n (2.1) 

n n 

Then 


*x(t) - *x°(t) 
m m 


< „0(e)v t e [0, T] 
p> 


m+1 


m+1 o 


m+1. 


x(t) - x u (t) < _0(e) V t e (0, T) 


n 


n 


m+1 , - m+1 o , . m+1 b,_ N 

x(t) - x (t) - x (t) 
n n n 


< “^(OV t[0, T] 


( 2 . 2 ) 

(2.3) 

(2.4) 


where x(t) and x°(t) are the solutions of (1.2) and (1.4), respectively, 

0(e) is a vector, each of whose components is of order e, T is an arbitrar- 

nH*l b 

ily large but finite constant, and ^x (t) is the boundary- layer correction 
derived in the following section. 

Equations (2.2) and (2.4) express uniform convergence and (2.3) expresses 
nonuniform convergence. Equation (2.1) ensures asymptotic stability of the 
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boundary layers for the simple linear system we are considering. (For the 
nonlinear case, an additional hypothesis is needed to guarantee that the ini- 
tial conditions are in the domain of influence of the stable equilibrium 
point. ) 

From (2.2) through (2.4) we note that if we are Interested only in the 

behavior of *x(t), or in the behavior of ^^(t) at large t (as compared 
m n 

with e), or both, then x°(t) will be a good approximation; it is only if 

1 

the behavior of x(t) for small t is of interest that the boundary- lay er 
n 

corrections "^x^t) need be computed. We also note that in the special case 
n 

that is diagonal, then (2.1) requires the diagonal entries to be 

n n 

negative. 

The theorem may be generalized in several ways. For examp.le, it holds 
if the elements of E are defined by 


e i j “0; i * j 


e ii = 1; i * « 


lim e i+l,l+l (c) 
£-*■0 e ^ U ) 


0; i > m 


(2.5) 


The special form of E as given by (1.3) is used solely for convenience. The 
results may also be generalized to the case where A is a function of e and 
t. In this case, we assume that A(e, t) has an asymptotic power series ex- 
pansion In e and write: 

A(e, t) - Ag (t) + Aj (t)e + A2(t)e 2 + ••• 

where the A^(t) are assumed continuous. In fact, because only the zeroth 
order terms of the expansions are computed, the results obtained here apply 
directly to the case where A is time- invariant but depends on e ; we need 
only replace A by A q in the final formulas. The t^me- invariant case was 
adopted only to simplify the algebra. 
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Although the theorem may be derived as a special case of more general 

results [c.f. 5], it is instructive and perhaps easier to pursue a direct 

proof of the theorem. In the following section we will prove the theorem and 

obtain an explicit representation for x (t). The main task is to show 

n 

that the boundary layers are asymptotically stable. As mentioned earlier, 
the analysis will employ the MAE method. 


III. ANALYSIS 

First consider the reduced system (1.4). Strictly speaking, this system 
(also called the outer system in the MAE method) is not required to satisfy 
ary boundary conditions but rather has its constants of integration determined 
by matching with the first boundary layer. However, for the zeroth order, 
this matching simply states that the constants of integration are to be deter- 
mined from *x°(0) = ^x [c..f. 10]. By (2.1), is nonsingular so that 

m m o 1 J J n n 6 

the solution of (1.4) may be written 
1 o 

x * exp 
m 


V - ‘a ^ 1 

mm m r. 


fm+l^m+1 -1 flH-l^l 
n n n m 


1 

x 

m o 


m+1 o 
x 
n 


m+1 A m+1 
A 

n n 


* m+1 1 
A 

n m 



‘a 1 

m m 


1 .m+1 
A 

m n 


fm+1 m+1 

-1 

raH 1 


[ n n 


n m 

l 


1 

x 

ill 0 


(3.1) 


In particular. 


1 

m 


x°(0) 


1 

x 

m o 


m+1 

n 


x°(0) 


m+1 ,m+l] 
A 

n n J 


m+1 1 1 
A x 
n m m o 


(3.2) 


Now consider the first boundary layer. The zeroth order term in the 
solution to this system is obtained by transforming the independent variable 
to n * j in (1.2) and taking the limit e -*■ 0. Denoting the dependent vari 
able by x*(t), this function satisfies 
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(3.3) 


. 1 1 
d x 
m 


dt 


1 


l 0 ; VO)) - l x 

m m mo 


. mfl 1 

m+1 m+1.1 11. m+1. m+1 nH-1 1 m+1 .m+2 m+2 1 uH-1 1, A , 

,A_ _x + _.,A _.,x + _.,A_ x ; _.,x (0) 


dT l 

m+2 


m+1 m m 


m+1 m+1 m+1 


m+1 n 


n 


m+1 


p+1 

m+l X o 


. m+2.1 11. m+2 .m+1 m+1 1 . m+2 .m+2 m+2 1 

0 * A x + A ,x + A x 

n n m m n m+1 m+1 n n n 


ml 2 m | 2 

Since *A^ is nonsingular, the solution of this system is 

1 1 1 
x * x 
m mo 


m+1 1 
m+l X 


m+2 1 
x 
n 


■(: 


m+1 m+l„l 1 

m+l X o m+1 m m X o 


exp 


fm+l^m+1 
n n 


m+2^m+2 1 

n n 


m+2 A m+2l~ 1 
n n 


("1 


.1 m+2 .m+1 m+lj.1 

m " n A m+l m+1 m 


m+lj.1 1 

m+1 m m X o 


1 

x 

m o 


. m+2. m+1 

+ nVl 


m+1 m+1 1 1 

m+l X o m+1 m m X o 


exp 


m+l^m+1 
n n 


m+2^m+2 1 

n n 


(3.4) 


where (1.1) was used and K is defined as 


V 

JJ-1 


JJ-1 


j A j+1 J+1 A j+1 
J n [ nn 


-1 


J +1 A J 


n J-lj 


J +1 D j+1 

n n 

J D J 

n n 


(3.5) 


By (2.1), we know that the solution of the first boundary- layer equation 
x^ is asymptotically stable. The matching condition of the MAE method re- 
quires that the solution of the first boundary- layer equation approaches, as 
-*• », the solution of the outer system evaluated at t ■ 0. The value of 
x* as -*■ *> may be found from the equilibrium point of (3.3) which is 

11,, 1 

m mo 


m+1 1 , , 
x (•) - 
n 


m+1 .m+1 
A 
n n 


-1 m+1 1 1 
A x 
n m m o 


(3.6) 
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Comparing (3.2) with (3.6) we see that the matching condition is satisfied, 
i.e. , 

x°(0) - x 1 ('*) 


and thus the reduced solution evaluated at the boundary is an asymptotically 
stable equilibrium point of the first boundary- layer equation, as required. 

Next, consider the i— boundary layer 1 < i < n - m. The zeroth order 
term of the solution to the i— boundary- layer system is obtained by trans- 
forming 



e 


*nd setting e « 0. Denoting the solution of the system by x*, 


(3.8) 


. 1 i 

W 

d t , 


1„ 1 i/AN 

m+i-1 0 ’ mfi-l x 


1 

mfi-l X o 


d ^i 1 

mfi mfi.l 1 i . mfi. mfi mfi i . mfi.mfifl mfifl i mfi mfi 

— 1 ” ..A . , . , ,x f ..A .. ,.x f ..A x; , ,x (0) ■ ,.x 

d t . mfi mfi- 1 mfi-1 mfi mfi mfi mfi n n mfi mfi o 


mf if 1 mfifl. 1 1 i , mfifl .mfi mfi i mfifl .mfifl mfifl i 

0 ” A . . . . ,x f A . .x f A x (3.9) 

n n mfi-1 mfi-1 n mfi mfi n n n 


nH"i i. 

We see that, in effect, the variables that are slower than .tx , the components 

m+1 

of j_jx*, have not as yet had time to move away from their initial values 

while those that are faster than **t*x*, the components of have al- 

mf I n 

ready left the region of influence of their initial conditions. Proceeding as 
before, the solution to this system is 
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We have now shown that all boundary layers are asymptotically stable and 
that (3.13) holds for all 1*0, l, ••*, n - m - 1. Application of elementary 
singular perturbation theory then completes the proof of the theorem (c.f. 1-4 
and 10] . It remains only to find a representation for the boundary- layer cor- 

_ . m+1 b, v 

rection terms x (t). 

n 

In the MAE method, uniformly convergent approximations for the fast vari- 
ables are obtained by forming composite solutions [c.f. 6, 7, and 10]. Since 
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fil+1 b 

there are many ways of doing thia, tha function *x (t) is not unique, In 
tha praaant analysis, wa adopt tha additive composition. To obtain tha addi- 
tive composition, the reduced and boundan -layer solutions are added and the 
common parts (l.e., tha terms that cancel out In the matching) are subtracted 
out. Tnls gives 

x“(t) - x°(t) + l x A (t) - n l 1 x i (0) (3.14) 

i-1 i-0 

From (3.2), (3.4), and (3.10), 

i x±(0) " iv 1 “ °» ■■■• n - m 

x (t) - x ; i - 1, ••*, n - m 
m mo 

Thus, for the slow variables 

*x a (t) - *x 0 (t) (3.15) 

m m 


and, as expected, no boundary- layer corrections are needed. 
For the fast variables, (3.2), (3.4), and (3.10) give 

jx i (0) - jx^; j s m + i 

1 i, . ^ 

jx (t) "j x o ; J < m + i 

Thus, from (3.14) 


jx a (t) 


j x0(t> 


1-m j , j-m-1 . . 

I l (t) - l ^x i ( 0 ) ; j 

i-1 3 i-0 3 


m + 1, 


(3.16) 


and the boundary- layer corrections are 


}x\t) 


l jx ( t ) — l jx (0); J • a + 1, 

i-1 J i-0 3 


n 


* 
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or, using (3.13), 


v<t> - 'i 


j-m 

l 

i-1 


jx^t) - 




n 


(3.17) 


so Chat only the exponential terms of the boundary- layer solutions appear in 
r x (t). From (3.17) and (3.10) the boundary -layer corrections may be writ- 
ten as 


"♦‘xNt) - tH ' 1 V(t) 

n n M n n 


(3.18) 


jjfl jb+1 

where R Q n is an n - m * n ~ m lower triangular matrix whose components 
q^ are given by 


l U 


0; J > i 




nrt-j nrfj„l 1 . , 

mfj X o + nrfj K nrt-j-l nrt-j-l X o’ J n “ m “ 

n . 1 n.l 1 

x + A . ,x ; j ■ n - m 

no a n n-1 n-1 o J 
nn 


_i- j ! nrt-j-f 1 V nrt-j+l 
q.. - - R J < J A J 

ij U n n 


1-1 


nrt-j +1 nrt-j 
n nrt-j 


nrt-j . nrt-j,, 1 1 ii . 

_ , .x + .^K .. . .. .x ); i > J , 

nrt-j o nrt-j nrt-j -1 nrt-j- 1 oj 1 


(3.19) 


where R*{A) denotes the i— row of A and where ^^(t) Is a vector whose 

n 

elements are scalar, stable, exponential functions, namely 


nrt-1 


V t) 


exp 


nrt-l^nrt-1 
ii n t 


exo 


nrt-2^nrt-2 c 
n n 


itrt-i_ nrt-J 

Ii 

n n 


nrt-i+1 nrt-i+1 i 
D e 

n n 


exp 


nn n-o 

e 


(3.20) 


where (3.8) was used. 
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IV. DISCUSSION AND APPLICATIONS 


In the previous section, we have obtained an approximate solution to the 

th th 

n — order system (1.2). This solution requires solving an m — order (m < n) 

system and a sequence of n - m scalar boundary -layer equations. Solution of 

the boundary- layer equations involves finding the determinants and inverses of 
m | ^ ^ 

the matrices X A , i = 1, * * * , n - m. Thus, the method in effect splits 
n n 

up a higher-order problem into several problems of lower order. By the 
theorem of section 2, the method gives a uniform approximation. 

Although the system (1.2) has a rather special form which may not be ex- 
pected to occur frequently in applications, the method developed here may be 
used to obtain approximate solutions to general linear systems. Suppose we 
have a system 

^ - Ay; y(0) = y Q (4.1) 


with (2.1) satisfied for some m < n, and, in addition. 


a 

— 

V 

> 

n-l D n-l 

> . . . > 

nri-l^mfl 

nn 


n n 


n n 


n n 


with the first boundary- layer solution decaying faster than any mode of the 
other solution. Then we proceed by premultiplying the left-hand side of 
(4.1) by E, thereby obtaining a problem of the form (1.2). We have thus 
inserted the ! 'small ,, parameter to create artificially a singular perturbation 
problem. Such a procedure has been adopted by several researchers in flight 
vehicle trajectory optimization and has proved to be quite successful [10, 11, 
12, 13, 14, 15]. In particular, Calise [11, 12] has utilized the multiple 
boundary- layer approach we are using here. The final step of the procedure 
is to set 6=1, its proper value. From (3.20) this gives 


• ■ i 

. L V < 
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nri-1 

n 


V(t) 


exp 


aH-lpiB+l 

n n 


exp 


nri-2j.nrt-2 
^ n n 


m+i^m+i 
n n 


mfi+l^m+i+l 
n n 


exp (a t) 
nn 


(4.3) 


The relations (4.2) ensure that the corrections to the reduced solution from 
the i— boundary layer decay faster than do those from the i - 1— 
boundary layer. 

In specific applications, experience with the system under consideration 
often shows the proper choice of slow variables and ordering of fast variables. 
In genet i* ■ the following procedure may be adopted. First, search the diagonal 
elements a^ of A for the one with the smallest value (i.e., the one with 


the largest absolute value among those that are negative) ; relabel the equa- 

for 


tions such that this element appears as a^. Next, compute 
all i and relabel the equations such that the a for which 


a , , a. 
ii in 

a . a 
ni nn 


a. . a. 
ii in 

a . a 
ni nn 


is largest appears as a n _j n _i* Then, continue on in this way 


until a .. .. is selected. After premultiplying by E, the results of 

m+i t nH"i 

section 3 may be applied. 


V. EXAMPLE 


Since a third-order system is the lowest order system that can exhibit 
all the features of the method, we will illustrate the method by using it to 
obtain an approximate solution of 
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“ a lll/l + a 12j/2 + « 1 3^/3 5 l/l (°) " a 

* a 2 12/ 1 + a 22j/2 + a 23l/3; i/2 (°) * 6 

- « 31 yi + a 32 t/ 2 + a 33i/3; 2/3 (0) ■ Y 


It is assumed that ( 2 . 1 ) and (A. 2 ) hold for m = 1 , i.e.. 


and, in addition, that 


3D3 ■ a 22 a 33 ” a 23 a 32 > 0 
3D3 - a 33 < 0 


I a 3 3 I > l a 22 a 33 ” a 23 a 32 


i a 22 a 33 ~ a 2 3 a 32 1 > l D 


where D is the determinant of A. The appropriate singular perturbation 
formulation of the problem is then 


= a^xi + a j 2 £ 2 + a i 3*^3 » ®i(0) = 


6 (H? - " a 2 l^l + a 22*2 + a 23^3J X 2 (0) = 
f. 2 = a 31^1 + a 3 2*^2 + a 33 ar 3» *3 (0) = 


We will solve this problem both by direct application of the MAE method and by 
employing the formulas derived in section 3 . 

First, the MAE method will be illustrated. The reduced problem is 

dx? 0 , 0 , 0 

= a ll^l + a l 2* c 2 + a 13^3 

0 = a 2 \Xi + a 22 x 2 + a 23^3 

0 = + a 32 :r 2 + a 33 :E 3 


with solution 
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l M ll 


Xi - C° exp 7 ~ t 
l M ll 

Xi ■ - C° exp 

Mi i 

0 Mi 3 n D 

*3 - HTT c “ p lsi 7 c i 


where is the minor of and C° is a constant to be determined by 

the matching condition. At t • 0, 

x9(0) - C° ; xS (0) - - Jr 2 - C°; «2(0) - ^ C° 

The first boundary- layer system is, to zeroth order, 

dxi _ 

dti 
.1 


dX2 1 . l . 1 

■ a 2 i*i + a 2 2*2 + a 2 3® 3 


0 = a 3 ]x| + a 32 X 2 + 333 X 3 


where t 1 * tie. The solution in terms of constants of integration 
c| * x}(0) and C 2 * X2<0) is 


_1 _ r l 
*1 C 1 




l 

*3 


M, 


rl + “L 2 rl 
L 2 + M n C 1 


f M ll Mi 2 , 

exp x 1 - 7 — C , 1 
l a 33 l J M 11 1 


a 32 

Ll A Mi 2 1 ] 


'Mn t ’ 

a 3 3 

[ C2+ m„ c ‘J 

| exp 

a 33 Tl . 


Mn 1 
+ XTT C 1 


At Ti * 0, 


a 3 2 
a 33 


1 Mi 2 1 ) M 13 1 

C2 + mtt Ci J + sir cl 


xj(0) * c}; X2 ( 0 ) « C2; x 3 (0) 

and for ti -*• », 

x[ (®) » c} ; X2 (°°) - - 
The zeroth order of the second boundary -layer system is given by 


M 12 r l. 1, v . M 13 1 

m m Cl * mTT 1 
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I 


Act 

dr 2 

.2 

d~2 

dt2 


« 0; 


0; 


x 2 (0) » a 


x^CO) - 8 


2 2 2 2 
■ a 31 x l + a 32*2 + a 3 3**- 3 > ^3(0) - Y 

where T 2 ■ t/e 2 . Note that since this is the last boundary layer it is re- 
quired to satisfy the initial conditions. The solution is 

xf ■ a 

x\ ■ 8 


* 3 


Y + ( a 31 a + a 32&) 


exp(a 33 t 2 ) - (a 31 a + a 32 B) 


As To -*■<*>: 


Xj (°°) = a; xfC 00 ) ■ 6; x 3 («) 


a 3 3 


(a 3 ia + a 3 2$) 


Next, the matching conditions are applied to determine the constants of 
integration. For matching between the first and second boundary layers we 
require that x^(°°) * x*(0) for i = 1, 2, 3; this gives Cj * a and C;? ■ 8. 
Similarly, the matching condition between the outer solution and the first 
boundary layer is x*(°°) = x°(0) for i * 1, 2, 3 which gives C° = a. From 
(3.15) and (3.17) the additive composite solution for y is 


y 1 
y 2 

y 3 


a exp 


H 


D 


12 

Mn 


Mu 


a exp 


t 


(Mu 


1 Mi 2 1 

+ 13 + Mrr a 


Mn 

exp — — t 
l a 3 3 


M 1 3 / D „ 

- — a exp It: — t 


'1 1 


Mil J a 3 3 


332 's 


Mil J 


exp 


(^■1 


Y + — ( a Sl a + a 32 5) 


exp (a 3 3 t) 


which provides a uniform approximation to the exact solution of the system. 
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We now solve the example by employing the results of section 3. The 
first step is to solve the reduced system, which we have done previously. It 
remains to obtain the boundary-layer correction terms. From (3.19), 


912 ■ 0 


2 2 11 
9ll " 2 X 0 + 2^1 i*o 


3 . 1 3.1 

922 " 3 X 0 + 3 A 2 

2 X 0 

i ffs a) -1 
921 “ ‘ R 1 |[3 A 3 

3 A 2^2 x 0 + 2 K 1 1 x 0 | 


From (3.5), 


2*1 



v 

2 A 3 


V 

3 A 3 


-1 


3.1 

3 A 1 




a 2 3 T7 7 a 31 
a 33 


±31 

Mil 


Ml 2 

Mil 


Thus, 


qi2 - 0 


9 ii 


6 + 


12 
1 1 


a 


922 


Y + 


1 

a 3 3 


( a 3 l a 


+ a22&) 


921 


- R 1 



‘32 


6 + 


Ml2 

Mil 



a 32 


a 33 


6 + 


Mu 


' 

a 

- 


From (3.20), 


§V(t) 


exp 


Mu t 

a 33 €. 


exp 


*33 77 
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Then from (3.18) 


n K 
3* <t) 


(■*fc -l-feij 


L 132. 


$ + 


H U 


a 3 3 i Mil J 


exp 


'Mll t 
a 3 3 


When 3 X^(t) is added to 3 X°(t) and 
tion results. 


+ 


e 


[ T + ik (a3i “ 

is set to 1, 


+ a 32fJ) 

. 

the desired approxima- 
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